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ABSTRACT 

We present 2.12-2.23 /im high contrast integral field spectroscopy of the extrasolar planet HR 8799 b. 
Our observations were obtained with OSIRIS on the Keck II telescope and sample the 2.2 ^tm CH4 
feature, which is useful for spectral classification and as a temperature diagnostic for ultracool objects. 
The spectrum of HR 8799 b is relatively featureless, with little or no methane absorption, and does not 
exhibit the strong CH4 seen in T dwarfs of similar absolute magnitudes. The spectrum is consistent 
with field objects from early-L to T4 (3 a confidence level), with a best-fitting type of T2. A similar 
analysis of the published 1-4 ^m photometry shows the infrared SED matches L5-L8 field dwarfs, 
especially the reddest known objects which are believed to be young and/or very dusty. Overall, we 
find that HR 8799 b has a spectral type consistent with L5-T2, although its SED is atypical compared 
to most field objects. We fit the 2.2 fim spectrum and the infrared SED using the Hubeny & Burrows, 
Burrows et al., and Ames-Dusty model atmosphere grids, which incorporate nonequilibrium chemistry, 
non-solar metallicities, and clear and cloudy variants. No models agree with all of the data, but those 
with intermediate clouds produce significantly better fits. The largest discrepancy occurs in the J- 
band, which is highly suppressed in HR 8799 b. Models with high eddy diffusion coefficients and high 
metallicities are somewhat preferred over those with equilibrium chemistry and solar mctallicity. The 
best-fitting effective temperatures range from 1300-1700 K with radii between ^0.3-0.5 i?j up - These 
values are inconsistent with evolutionary model-derived values of 800-900 K and 1.1-1.3 i?j up based 
on the luminosity of HR 8799 b and the age of HR 8799, a discrepancy that probably results from 
imperfect atmospheric models or the limited range of physical parameters covered by the models. The 
low temperature inferred from evolutionary models indicates that HR 8799 b is ^400 K cooler than 
field L/T transition objects, providing further evidence that the L/T transition is gravity-dependent. 
With an unusually dusty photosphere, an exceptionally low luminosity for its spectral type, and hints 
of extreme secondary physical parameters, HR 8799 b appears to be unlike any class of field brown 
dwarf currently known. 

Subject headings: planetary systems — stars: individual (HR 8799) — techniques: image processing 



1. INTRODUCTION 

High contrast adaptive optics (AO) observations on 
large telescopes are beginning to explore the outer ar- 
chitecture ( ^,10 AU) of plane t ary systems for the first 
time (e.g . iBiller et al.l [20071: iLafreniere et al.l l2007at 
iLiu et ail 120091: iChauvin et al.l I2010D . The recent 
discoveries of giant planet system s around HR 879 9 



dMarois etaU l2008bD. Fomalh a ut dKalas et~ 
P Pic (jLagrange et al.l 120091 : Lagrange et al 



and 1RXS J160929.1-210524 ( Lafrcni ere et al J 1200a 



2008) 



20101) 



ILafreniere et al.ll201Q[) by direct imaging have shown that 
planets exist at large separations (^10-330 AU) and are 
within reach of current telescope capabilities. The ad- 
vantages of direct imaging compared to the radial veloc- 
ity technique are threefold: it is the only way to study 
planets on wide orbits; it allows for estimates of the true 
planet mass through luminosity measurements (rather 
than minimum masses as is the case for non-transiting 
radial velocity planets), and it enables follow-up photo- 
metric and spectroscopic observations to study extrasolar 
planetary atmospheres. 
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Multiwavclcngth observations of extrasolar planets can 
be used to infer their physical and chemical proper- 
ties and test our understanding of planetary atmo- 
spheres through direct comparisons with theoretical at- 
mospheric models. Transmission spectroscopy and sec- 
ondary eclipse observations of transiting cxoplancts have 
led to the detections of several molecular species in 
the atmospheres of h ot Jupiters, namely H 2 Q, CO2, 
CO and CH 4 (e . g.. iTinetti et al.l l2007t iSwain et al.l 
120081 : ISwain et all 12009 1 These short-period transit- 
ing gas giants are also providing valuable tests of irra- 
diated planetary atmospheric models a nd have shown 
evidence for temperature inversions (e.g.. iKnutson et al.l 
l2009t iSpiegel et ail 2009f) , global atmospheric dynam- 
ics ((Showman et al.ll2009f). and non-equilibrium carbon 
chem istry (e.g., iCooper fc Showman! 120061 : ISwain et al.1 

l2oTnh . 

At much larger separations of tens to hundreds of 
AU, observations of young (<100 Myr) substcllar com- 
panions near and below the brown dwarf /planetary- 
mass limit (~13 Mj up ) are revealing the physical prop- 
erties of very low-mass objects throug h direct imag- 
ing and spectroscopy (e.g ., AB Pic B: IChauvin et al 
20051 iBonnefov et alJl2010t GQ Lup B:lNeuhauser et al 
20051 iMcElwain et al.l 120071 "~ 
ble 1 of iZuckerman fc Songj 



Lavigne et al. f l200¥ see Ta- 
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For example, atmo- 
spheric model comparisons to near-infrared spectroscopy 
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of the planetary-mass companion 2MASS J1207-3932 b 
(|Chauvin et al.l [20041 ) is demonstrating the important 
role th at dust may play in the atmospheres of you ng gas 
giants (jMohantv et alJl2007t [Patience et al.ll2010t ). 

The system of three directly imaged giant planets or- 
biting HR 8799 presents an excellent opportunity to 
investigate the properties of extrasolar planetary at- 
mospheres at young ages. HR 8 799 (= HP 218396 , 
V342 Peg) is a youn^l (60±^° Myr; IMarois et~aHl2008bl ), 
intermediate- mass (1.5 ± 0.3 M Q ; IMarois et al.ll2008bl) . 
A-typ^3 star with an Hipparcos par alla ctic distance 
of 39.4 ± 1.1 pc (Ivan Leeuwenl 120071 ). IMarois et all 
(|2008bD discovered three faint common-proper motion 
companions with projected separations of 68 AU ("b"), 
38 AU ("c"), and 24 AU ("d") and evolutionary 
model-derived masses of 7+_\ M Jup , 10±| Mjup, and 
lO^g Mj up , inferred from their luminosities and the age 
estimate for HR 8799. Subsequent an alysis of previ- 
ous i maging data from Ke ck II in 2007 (iMetchev et al.1 
I2009T) . Subaru in 2002 (|Fukagawa et all I2009D. and 
the H ubble Space Telescope in 1998 ( Lafreniere et all 
120091) revealed several pre-discoveries of the HR 8799 
planets and confirmed their common proper and or- 
bital motions. Dynamical modeling of this sys- 
tem suggests a small (but non-zero) inclination and 
two- or three-body resonances a mong the planets 
(iGozdziewski fe Mi gaszcws kH 120091 i Rcidcm cister et all 
120091: iFabrvckv fc Murrav-Clavll2010bir Thermal emis- 
sion has been detected from a warm inner debris disk 
around HR 8799 at ~10 AU and an outer massive disk 
at >100 AU, implying an architecture analogous to our 
Solar System with giant planets surrounded by "exo- 
aster oi d" and "exo-Kuiper" belts dSadakane fe Nishidal 
l98l iZuckerman fe Sonej 12001: iMoor et all 120061: 



Williams fe Andrews|l200alRhee et ai.ll2007t iChen et all 
2009HSu et all 120091 : iReidemeister et alJl2009f ). ~ 



The atmospheric properties of the HR 8799 plan- 
ets have been investigated in several studies. All 
three planets exhibit redder near-infrared colors than 
field T dwarfs with similar absolute magnitudes, which 
suggests th e presence of dusty clouds in their at- 
mospheres (jMarois et all I2008bl : iMetchev et"all 120091 : 



3 Recently, Moya et al. (2010a) inferred an age of up to sev- 
eral Gyrs for HR 8799 using asteroscismology, which would imply 
that the HR 8799 companions are brown dwarfs (^30-40 A/j U p 
at 1 Gyr; [Baraffc et al. 20Q3). However, the Frequency Ratio 
Method l|Mova et al.l 120051 ) used in that study is on ly accurate 
for a range rotational velocities (Suarcz et al. 2005) , which for 
HR 87 9 9 limits the applicability t o inclinations >36° (Mova et al. 
I2010bl ). Rcidcmcistcr ct al. (2009;) infer an inclination between 20- 
30° based on a stability analysis of the companions and combin- 
ing the range of true rot a tional velocities with vs'mi. Moreover, 
Fabrvckv & Murrav-Clay (2010a) find that the orbiting compan- 
ions must have masses <20 Mj up based on a stability analysis. Al- 
together this points to a younger age for the system. Unfortunately, 
searches for wide stellar companions to HR 8799 which could fur- 
ther hel p constrain the age have not produced any promising can- 
didates flClose fc Males 20091: Iffinz et al.llMiOl) . Addi tionally, the 
subsola r internal metallicity of HR 8799 inferred by IMova et al.l 
l|2010bl ) may not be accurate because of the most probable incli- 
nation of the system. 

4 HR 8799 is a 7 Doradus variable star and belongs to t he class 
of chemically peculiar A Bootis stars. IGrav fc~K avc (1999) infer a 
spectral type of kA5 hFO mA5 V A Boo based on detailed spec- 
troscopic analysis, where "k", "h", and "m" refer to the Ca II K, 
hydrogen, and metallic line strengths and profiles, respectively. 



lLafreniere et all 120091 ). Additionally, non-equilibrium 
CO/CH4 chemistry in t he atmospheres of HR 8799 c and 
d has been suggested bv lHinz et all (|2010l) based on mid- 
infrared (3.3 (im, L'-band, and Af-band) photometry. 
Likewise, iJanson et all (|2010T l presented a low signal-to- 
noise (S/N) 3.88-4.10 /im spectrum of HR 8799 c and 
suggested that the disagreement with atmospheric mod- 
els based on chemical equilibrium may hint at the influ- 
ence of non-equilibrium chemistry. 

Methane forms at temperatures <1400 K (e.g., 
iFortnev et al.ll2008l ) , so under equilibrium conditions the 
spectroscopic signature of CH4 provides an independent 
diagnostic of effective temperature for ultracool objects. 
The presence of CH4 absorption in the near-infrared is 
also the defining spectral signature of T dwarfs and there- 
fore is a useful tool for spectral classification. Medium- 
band photometry on and off the 1.6 fim methane band 
with the CH4S and CH4I filters shows no evidence for 
methane absorption in HR 87 99 b and c, althoug h it 
may be present in HR 8799 d (jMarois et alll2008bh . In 
this work we focus on the 2.2 /xm absorption band of 
HR 8799 b, which is also sensitive to non-equilibrium 
CO/CH4 chemistry at low temperatures and low surfac e 
gravities (|Hubenv fc Burrowsll2007t IFortnev et alll2008t l . 

We describe our adaptive optics (AO) observations 
with OSIRIS at Keck II and our data reduction proce- 
dure in Sj2]and ^3j In $4] we compare field L and T dwarfs 
to our spectrum and previously published photometry for 
HR 8799 b, and we do the same with atmospheric models 
in fJSJ Finally, we discuss our results in SjB]and we provide 
a summary in ^7] 

2. OBSERVATIONS 

High contrast AO imaging observations are limited 
by quasi-static speckle noise at small angular sepa- 
rations (<l-2"), which results from slight imperfec- 
tions in telescope o p tics and AO wavefr o nt correction 
dRacine et all 119991: IMarois et all 120001: IMarois et all 
I2008at lOppenheimer fc Hinklevl l2009ll These speckle 
patterns pose many problems for directly imaging faint 
targets near bright stars. Speckle noise is correlated and 
adds coherently, so longer exposure times do not im- 
prove contrast performance in this regime. The speckles 
themselves can resemble real astronomical objects, mak- 
ing it difficult or impossible to distinguish between noise 
and true signal in a dataset. Additionally, speckles are 
time-depen dent phenomena spanning a range of dura - 
tions (e.g., iMacintosh et all 120051 : iHinklev et all 120071 ). 
Because of the temporal evolution of speckle patterns, 
no single observation of the point spread function (PSF) 
at a given time can accurately represent the PSF for ob- 
servations taken over a long period of time. 

As a result of these difficulties, there has been a con- 
certed effort to develop observing and reduction tech- 
niques to overcome the limitations imposed by speckle 
noise. One ver sion of the spectral differential imaging 
technique fSDl lSparks fc Fordll2002l: iThatte et al.ll20071) 
takes advantage of the wavelength-dependent nature of 
speckle positions to distinguish speckles from bona fide 
objects, which remain stationary across all wavelengthsQ 

5 This version of SDI makes the fewest assumptio ns about the 
planet's spectrum. An earlier, simpler version of SDI IISmitblll987k 
IRacine et all 119991 : IMarois et aT] 2000; Biller ct al. 20071) involves 
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In principle, integral field spectro graph (IFS) observa- 
tions can naturally utilize SDI (iSparks fc Fordl l2002t 
Berton et all 120061 iThatte et alj 120071: lAntichi et al.l 
20091) . Yet to preserve the flux from real objects, SDI 
requires a fractional change in wavelength that is large 
enough to avoid self subtraction after de-magnification. 
SDI with IFSs is therefore more practical for observations 
with broad wavelength coverage (i.e., broad-band filters 
rather than narrow-band filters) . Here we make use of a 
variant of the angular differen tial imaging method (ADI; 
[Tiv][200l |M arois et all 120061 ). which normally relies on 
rotation of the field of view when the telescope rotator 
is turned off. This speckle subtraction technique makes 
it possible to distinguish between the speckle pattern, 
which remains stationary, and real objects, which rotate. 
Our implementation of ADI is different, however, as we 
leave the rotator on to ensure that HR 8799 b remains 
in the same position on our detector (to avoid the planet 
rotating out of the narrow field of view of our obser- 
vations; see below). In this version of ADI the planet 
remains stationary while the speckles rotate. 

Our observations of HR 8799 b were carried out 
with the OH-Sup pressing InfraRed Integral-field Spec- 
trograph (OSIRIS: lLarkin et al.ll2006l ) with natural guide 
star AO at Keck II on 2009 July 21 UT. The weather 
was clear with a seeing estimate of ~0f.'7 in the opti- 
cal. OSIRIS IFS observations result in data cubes with 
spectral information over a large number (between ~1- 
3xl0 3 ) of spatial positions. The longer-lived portion of 
the speckle pattern in the OSIRIS data cubes can be used 
to register the individual images in an ADI sequence. 
Ideally we would obtain broad wavelength coverage over 
a large field of view to maximize the number of speckles 
in our observations, but there is a tradeoff between spec- 
tral coverage and field of view with IFS instruments as a 
result of limited space on the detector. We chose the Kn3 
filter (2.121-2.229 /Ltnfi Figured]) with the smallest plate 
scale (0'.'02) to target the 2.2 /im CH4 band with a mod- 
erate field of view (~CK'96xl'/28) while simultaneously 
well-sampling the PSF. Our filter choice was also based 
on AO performance, which is best in the if-band, and 
the fact that better CH4 line lists exist at 2.2 fim than 
at 1.6 (im, thereby enabling more accurate constraints 
on physical parameters using atmospheric models. The 
spectral resolving power (R = A/AA) of the OSIRIS data 
varies as a function of lenslet geometry (which depends 
on the filter), wavelength, and spatial position. At the 
location of HR 8799 b in our observations, we estimate 
R ~ 43000 

We obtained 18 consecutive observations of HR 8799 b 
with exposure times of 300 s each for a total integra- 
tion time of 90 min. The airmass varied from 1.16 to 
1.00 throughout our observations, spanning 13° of the 
rotation of the PSF in our field of view. Immediately 
afterwards we obtained observations of the A0V stan- 
dard star HD 208108 for telluric correction. We took sky 

simultaneously imaging in two neighboring wavebands over which 
plant-star contrasts significantly vary, e.g. at the 1.6 (im CH4 
absorption band present in low-temperature objects. 

6 These wavelengths correspond to the half-power points for the 
Kn3 filter. The OSIRIS pipeline excludes wavelengths outside this 
region in the reduction. 

7 Based on data in the O SIRIS User's Manual, v. 2. 3 
< |http: / /irlab. astro. ucla.edu/osirls}. 
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Fig. 1. — Keck II/OSIRIS Kn3 filter transmission profile. Over- 
plotted are IRTF SpeX/prism spectra of L and T dwarfs show- 
ing the increasing strength of the 2.2 fim CH4 feature at cooler 
temperatures. The spectra were obtained from the SpeX Prism 
Spect ral Library and were originally published by Burgasscr et al. 
BOTH 2MASS J02271036-1624479: LI). IReid et alJ J2006l . 
2MASS J10101480-0406499: L 6). ILooper et al.l 1120071, SDSS 
J120747.17+024424.8: TO), and IBurgasser et al.l 1 120041. 2MASS 
J15031961+2525196: T5, 2MASSI J0415195-093506: T8). All data 
are normalized to 2.1 /im and offset by a constant. 

frames before and after the science and standard star 
sequence. The FWHM of the standard star is 2.4 spax- 
els at 2.2 ^m, or 0'.'048 with the 0'.'02 plate scale (i.e., 
diffraction- limited) . 

3. DATA REDUCTION 
3.1. LOCI Processing 

Initial data reduction was perfo rmed using the OSIRIS 
data reduction pipeline (v. 2. 2; iKrabbd 120041) . The 
OSIRIS pipeline corrects for detector artifacts, performs 
sky subtraction and flat fielding, wavelength-calibrates 
the data, and assembles the 2D spectra from each lenslet 
into 3D data cubes. Every resulting data cube contains 
3063 spectra, each with 433 spectral channels and a spa- 
tial geometry of ^48x64 "spaxels," or spatial pixels. An 
example of a median-combined data cube from our ob- 
servations is shown in Figure [2] (left panel). 

The limited field of view of our OSIRIS configuration 
made registering the individual data cubes more difficult 
than for standard high contrast imaging observations be- 
cause the OSIRIS field of view encompassed only a por- 
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Fig. 2. — Left: Example of a median-combined OSIRIS data cube 
from a single 5 minute observation of HR 8799 b. The planet is 
mostly indistinguishable from the speckle noise in this image. Con- 
secutive images show the speckles rotating but the planet fixed in 
the same location. The location of HR 8799 A is ~0.8" below the 
edge of the array. Right : Fully processed image of HR 8799 b. The 
planet (marked with an arrow) is clearly visible after the image re- 
duction using LOCI. This image is the median-combination of 18 
individual observations for a total integration time of 90 minutes. 
LOCI has been applied in two annular subsections using the pa- 
rameter s from j|3l The imag es are displayed with an asinh intensity 
stretch IILuptoriet al.H2004r ). Note that the object ~0'.'6 below the 
planet is an artifact from poor speckle subtraction. 

tion of the stellar PSF and did not include the peak. 
(HR 8799 was placed ~0'/8 below the edge of the array.) 
Slight drifting caused the image positions to move on 
a 1-2 spaxel scale over the course of our observations. 
To register the data cubes we chose an observation near 
the middle of the ADI sequence and then used the rela- 
tive rotation angle of the remaining observations to in- 
dividually derive the x and y pivot coordinates about 
which the telescope rotated. The rotation angle between 
exposures was computed from the FITS header infor- 
mation and remained a fixed parameter in our analysis. 
We binned each data cube into 10 spectral channels and 
solved for a single pivot point position for each obser- 
vation. This was done in an iterativ e fashion using the 
AMOEBA algorithm (jPress et al.ll2007f ) by minimizing the 
rms of each de- rotated residual data cube (i.e., the fixed 
cube in the middle of our ADI sequence minus the cube 
being considered). The rms was computed in an over- 
lapping rectangular region of the residual data cube and 
included several bright speckles. A visual inspection of 
each de-rotated data set confirmed the proper alignment 
of bright speckles^ We tested the registration with sev- 
eral choices for the number of spectral bins and found 

8 Although the speckles were well aligned using our registration 
technique, there was a large variation in the inferred pivot point 
positions (usually differing by 5-20 spaxels, but even larger for a 
few of the exposures). The reason for this discrepancy is that the 
pivot point position is highly sensitive to small adjustments when 
aligning two exposures using this technique. The impact on the 
extracted spectrum was minor when the pivot point and rotation 
angle were fixed in the reduction; the main effect was to increase 
the spectral measurement errors, but the shape of the spectrum 
was preserved. 



TABLE 1 
HR 8799 b Spectrum (OSIRIS/Kn3) 
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that 10 bins adequately compensated between S/N and 
number of spectral channels. 

We performed speckle subtraction on our de-rotated 
data cubes using the "locally optimized combination 
of images" (LOCI) algorithm (jLafreniere et all l2007bt 
Figure [2j right panel). For each image, LOCI solves 
for the optimal linear combination of all the other im- 
ages to build a reference frame for PSF subtraction. 
We median-combined each data cube into 10 spectrally 
binned frames and independently treated each binned 
spe ctral channel as images i n an ADI sequence. Follow- 
ing eniere et al.1 (|2007b| ). we experimented with the 
different LOCI optimization parameters that control the 
shape and size of the optimization annular subsection 
(g and Na), the size of the subtraction subsection (dr), 
and the minimum displacement distance from the science 
target to avoid self-subtraction (Ns). We systematically 
varied each LOCI parameter and the number of bins to 
study the influence on the resulting spectrum. Adjusting 
the values of N A (from 100 to 400), g (from 0.5 to 2.0), 
and dr (from 2 to 10) had little effect on the spectrum. 
Values of N 5 below 1.0 (in units of PSF FWHM) resulted 
in a smaller signal from the planet due to increased self- 
subtraction, but values of Ns between 1.0-2.0 had only 
a minor impact. Increasing the number of spectral bins 
resulted in a noisier spectrum, but the shape and am- 
plitude of the spectrum remained unchanged. We used 
Na=300, g=1.0, dr=10, and iV<s =1.0 for our final spec- 
trum. 

After speckle subtraction, we performed aperture pho- 
tometry on each individual data cube to extract the spec- 
tra of HR 8799 b. Although the sky background should 
be zero as a result of using LOCI, we tested aperture 
photometry with and without sky subtraction as a pre- 
caution. ADI+LOCI produces a slight negative trough 
on either side of the planet spectrum as a result of de- 
rotation and subtraction. These troughs may affect sky 
subtraction using sky values computed from an annulus 
so we also tried sky subtraction using a sky value com- 
puted from a nearby box; the shape and amplitude of 
the extracted spectrum were preserved in all cases and 
the final flux measurements were all in agreement. We 
chose an aperture radius of 2 spaxels with no sky sub- 
traction for our final extracted spectrum. The final spec- 
trum was assembled by scaling each of the 18 spectra to 
the median-combined spectrum and then computing the 
median of the scaled spectra. 

To telluric correct the spectrum of HR 8799 b we first 
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Fig. 3.— Telluric-corrected OSIRIS spectrum of HR 8799 b. The 
mean S/N per spectral bin is 11. The data are normalized to 
2.1265 fim and the error bars are derived from our analysis of in- 
jecting and extracting artificial planets (see [|3]and FigureHJ. Weak 
CH4 absorption is suggested by the data from the negative slope 
between 2.15-2.22 (an, but strong 2.2 fim absorption is absent. We 
note, however, that the apparent dip at 2.18 fim is offset from the 
nominal 2.20 fim CH4 bandhead. 

extracted the spectrum of the standard star by spectrally 
binning the standard star data cubes, performing aper- 
ture photometry with an aperture radius of 10 spaxels 
(the FWHM is ~2.4 spaxels), and combining the indi- 
vidual standard spectra in the same fashion as we did 
for the planet. We then used the xt ellcor_basic pack- 
age in Spextool (|Cushing et al.|[2004}) to divide the planet 
spectrum by the standard spectrum and multiply by a 
blackbody with a temperature of 10,000 K (correspond- 
ing to an AO star) to restore the continuum shape. 

Our telluric-corrected spectrum is displayed in Figure 
[3] and the measurements are listed in Table [T] Methane 
can influence the 2.15-2.23 fim spectral region, with a 
particularly promine nt absorption feature at 2.20 /im 
i|Cushing et aL|[2005l ). Strong 2.2 fim CH4 absorption 
is noticeably absent in our spectrum, although a weak 
methane feature may be present as a negative slope from 
2.15-2.22 fim and diminished flux from 2.20-2.22 fim. 

3.2. Determining Measurement Errors 

Many sources of random and systematic errors can con- 
tribute to the uncertainties in our OSIRIS spectrum, in- 
cluding speckle noise, the use of a fixed radius in aperture 
photometry, and self-subtra ction of the science t arget as 
a result of ADI+LOCI fe.g.. lLavigne et al.|[2009l) . To as- 
sess the systematic and random errors in our spectrum, 
we injected and extracted fake planets in our data. For 
the artificial planets we used the standard star scaled 
to the amplitude of HR 8799 b. We injected the fake 
planets in six locations in the original individual images, 
reduced the full data set, and extracted the planets in 
the same fashion as we did with the science spectrum. 
To minimize the impact on the construction of the ref- 
erence PSF, we input one fake planet into the data at 
a time. The ratio of the input to output spectra of the 
fake planets varied on the level of ~5-10% depending on 
the location of the injected planet (Figure @|. In some 
locations systematic trends as a function of wavelength 
were observed, especially in the region closer to HR 8799 



where the speckle subtraction was worse. Overall, how- 
ever, a linear fit to the mean and rms of all six ratio 
spectra indicate that no correction is warranted as the 
slope of the fit is statistically consistent with zero. 

When random noise is the dominant source of error, 
the best estimate of the actual flux at each wavelength is 
the mean value, and the best estimate of the uncertainty 
in that measurement is the standard error of the mean 
{o~m = cr/v / iV, where TV is the number of spectra). Here 
we use the median flux value at each wavelength instead 
of the mean to avoid outliers in the data. In this case the 
expected uncertainty from photon noise is the standard 
error of the median (which for large sample sizes and nor- 
mally distributed dat a is i/(7r/2) x aj\/~N ~ 1.25 om\ 
IHoio fc Pe arson 1931). However, our analysis of the fake 
planet spectra indicate that the standard error of the me- 
dian of the real planet spectra underestimates the actual 
measurement uncertainty. 

To compute the uncertainties from our fake planet sim- 
ulations, we first computed the spectra for each of the 
six fake planets by taking the mean spectrum of 18 mea- 
surements from the individual images. After the spec- 
tra of the 6 fake planets were scaled to a common level, 
the standard deviation of the six spectra at every wave- 
length was adopted as the true standard error of the 
mean. When we compared these uncertainties from the 
fake planets to the standard error of the median from 
the 18 real spectra of HR 8799 b (and accounting for the 
factor of 1.25 between the two datasets), we found that 
the uncertainty derived from the spectra of HR 8799 b 
was underestimated by a factor of 1.5-2.5. We therefore 
adopted the standard error of the mean from the fake 
planet analysis for our HR 8799 b spectroscopic uncer- 
tainties. The final S/N per spectral bin of our OSIRIS 
spectrum varies between 7-18 with a mean value of 11. 

4. COMPARISON TO FIELD BROWN DWARFS 

In many ways brown dwarfs are massive analogs of gi- 
ant planets, sharing similar radii, effective temperatures, 
atmospheric physics, and cooling histories. Brown dwarfs 
can therefore serve as natural reference objects to guide 
our understanding of giant planet properties. Here we 
compare both our OSIRIS spectrum and previously pub- 
lished photometry of HR 8799 b to near-infrared spectra 
of L and T dwarfs from the Infrared Telescope Facil- 
ity (IRTF) SpeX Prism Spectral Library^ The spectra 
have resolving powers between 75 and 120 and were re- 
duced and telluric corr ected using Spextool data reduc- 
tion p ackage for IRTF (jVacca et al.ll2003t ICushing etahl 
120041 ) . When available we use optical spectral types for 
L dwarfs, and for all other objects we use near-infrared 
types. Objects with peculiar near-infrared classifications 
are omitted unless they have normal L-type optical spec- 
tral types. Altogether we use 238 L and T dwarf spectra 
in our analysis, consisting of 92 objects with L-type op- 
tical classifications, 56 with L-type near-infrared classifi- 
cations, and 90 with T-type near-infrared classifications. 

4.1. 2.12-2.23 fim Spectroscopy 

We fit the L and T dwarf prism spectra to our OSIRIS 
spectrum by binning them onto the same wavelength grid 

9 Maintained by_ Adam Burgasser at 

http://www.browndwarfs.org/spexprism 



6 




2.18 

X (|im) 

Fig. 4. — Results of injecting and extracting six artificial planets 
in our data. The fake planets were constructed from the stan- 
dard star observations, scaled to the amplitude of HR 8799 b, in- 
serted into all the individual images, and reduced using the LOCI 
algorithm. The ratio of the input to output flux are shown at 
each spectral bin (gray lines). A weighted linear fit to the mean 
and rms of the spectra (filled circles with error bars) yield a best- 
fitting slope consistent with zero (red solid line), suggesting that 
no systematic correction to the extracted spectrum of HR 8799 b 
is necessary. The 1 a confidence interval for the linear fit is shown 
as red dashed curves, and the equation displays the best-fitting 
slope, offset, and parameter uncertainties. A reduction in the ex- 
trac ted flux (F ou t /Fj n ) is a characteristic of the LOCI algorithm 
(see ILafreniere et al.ll2007bT ) . 

as the OSIRIS data and scaling each prism spectrum k by 
a factor Cfc, which is found by minimizing the x 2 statistic 



Xk 



E 



(1) 



where fi and er/^ are the measured flux density and un- 
certainty of HR 8799 b and Tk,% is the prism spectrum 
flux density measurement at each wavelength i. Equating 
the derivative to zero and solving for the scaling factor 
we get 



(2) 



The best-fitting spectral type to our HR 8799 b 
spectrum is T2 (Figure [5]). Objects between L8 and 
T3 yield relatively small % 2 values (<20 for 9 d.o.f.). 
Those later than T4 produce poor fits because of their 
strong methane a bsorption. We use the A\ 2 method 
(jPress et al.ll2007l ) to determine which templates are sta- 
tistically consistent with our data. We consider models 
with x 2 values <35.6 as being in agreement with our 
spectrum, which corresponds to a A^ 2 value with an in- 
tegrated probability of 99.73% (3 a) for 9 degrees of free- 
dom (25.3) plus the minimum x 2 value of the fits (10.3). 
This corresponds to spectral types earlier than T4. For 
an integrated probability of 68.3% (A^ 2 for 9 d.o.f. = 
10.4), our OSIRIS data are consistent with spectral types 
of -L4-T3E3 

10 The best-fitting object results in a reduced x 2 value (x 2 /^; 
for i/=9 degrees of freedom) of 1.15 using the standard error of the 
mean uncertainties from the fake planet analysis. The reduced x 2 
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Fig. 5. — Results of fitting 238 field L and T dwarf spectra 
to our OSIRIS spectrum of HR 8799 b. Top: x 2 values as a 
function of spectral type. The best-fitting spectral type is T2, 
although objects with spectral types earlier than T4 are consis- 
tent with the spectrum of HR 8799 b (see Q. Objects below 
the dotted line (x 2 =35.6) are consistent with our data at the 
3 a level. Bottom: Best-fitting spectra of field L6, T2, and T5 
dwarfs compared to our OSIRIS spectrum. Fro m top to bottom 
the s pectra represent 2MASS J1036 530-344138 HBurgasser et~aT1 
|20T(T) , SDSS J075840.33+ 324723.4 iTB urgasscr et al. 2003), and 
2MASS J0755480+221218 l IBurgasser e t al. 2006j) with reduced x 2 
values (9 degrees of freedom) of 2.30, 1.26, and 7.23, respectively. 

4.2. 1.1-2.4 M 771 Photometry 

We perform a similar analysis using previou s ly pub - 
lished photometry of HR 8799 b. iMarois et all (j2008bT ) 
presented J, CH^s, H, CH4I, K$, and L' pho- 
tometry of HR 879 9 b in their discovery paper. 
ILafreniere et~aLl (12001 published HST F160W p h otom - 
etry. iFukagawa et all 1)20091) and iMetchev et all (|2009f ) 
-ff-band photometry and lHinz et al.l ( 2010D an V detec- 



valuc using the standard error of the median uncertainties from the 
18 scaled spectra of HR 8799 b result in a reduced x 2 value of 2.39, 
which is further evidence that adopting the standard error s fro m 
the extracted spectra underestimate the uncertainties (sec i|3.2l ). 
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ing iCushing et al.l (|2008[) we compute a G statistic for 
each prism spectrum k: 



Fig. 6. — Empirical comparison of field L and T dwarfs to 1-4 jxm photometry of HR 8799 b from IMarois et al.l | |2008bl ). Top Left: 
Goodncss-of-fit statistic G' (Equation[3]l as a function of spectral type for fits to the J, CH4S, H, and Kg bands. Photometry is synthesized 
from 238 objects in the SpeX Prism Spectral Library. Ultracool dwarfs with spectral types of L5-L8 provide the best fits. See Figure[7]for 
the best-fitting individual objects. Bottom Left: Best-fitting L6, T2, and T5 spectra of field dwarfs scaled to the photometry of HR 8799 b 
at 10 pc. Mid- to late-T dwarfs provide poor fits to the NI R SEP of HR 8799 b as a result of brighter J-band and fai nter i\-band flux . 
The spectra represent 2MASS J21481628+4 003593 (L6. red; [Topper et al.ll200Sft , SDSS J143553.25+112948.6 (T2, blue MChiu eTTaTEXKjft . 
and 2MASS J23312378-4718274 (T5, green; ITjurgasser et al.ll2004l) . Toy Right: Goodness-of-fit statistic G" (Equation Q3 a s a function of 
near-i nfrared spectral type for fits of 40 ultracool dwarfs with J-, H-, Kg-, and L'-band photometry from the compilation in lLeggett et all 
112010 . their Table 3). Mid- to late-L dwarfs fit the 1-4 /im photometry of HR 8799 b the best. Bottom Right: Best-fitting L6, T2, and T5 
objects. From top to bottom the spect ra represent 2MASS J08251968+2115521, SDSS J075840.33+324723.4, 2MASS J15031961+2525196 
(for details see references in Table 3 of Lcggctt et al. 2010). 

tion with upper limits at 3.3 /zm and M-band. We use 
only the Marois et al. J, CH4S, H, and Ks photometry 
in this analysis to avoid overweighting the iJ-band region 
of the spectrum in our fits. 

We synthesize photom etry (e.g., Equation 5 of 
iTokunaga fc Vaccal 120051) from each prism spectrum 
of the field dwarf sample using the J, H, and Ks 
bands from the Mauna Kea Observatory (MKO) filter 
conso rtium ([Simons fc Toku naaa 2002; Tokunaa a et al.l 
20021) and the zer o poin t Vega flux densities from 
Tokunaea fc Vaccal (120051) . For CH 4 s, we use the 
Keck II/NIRC2 filter transmission curve (J. Lykc, 
private communication) and the flux-calibrated Vega 
model spectrum from Spextool (jCushing et al.1 120041 ) 
to compute a zero point flux density for that filter 
(5.31xl0- n W/m 2 /^m). 

We use a modified form of the x 2 statistic to assess 
goodness-of-fits for our photometric comparison. Follow- 



G, 



i=l 



a 



(3) 



f,i 



where Wi is the weight applied to each photometric point 
i and (J~k,i) is the monochromatic flux density in each 
bandpass. The scaling factor is then 



Cu 



w i (F) 



(4) 



We experimented with equal weights (u>i=1.0) for each 
filter, in which case G k becomes the usual xt statistic, 
and weights proportional to the width of each spectral 
bandpass such that ^iUj=l, defined as 



8 



(5) 



where n is the number of filters used in the fit. This 
results in weights of 0.18, 0.33, 0.13, and 0.36 for the J, 
H, CH4S, and K$ bands. We use this weighting scheme 
in our analysis, although using equal weights produces 
similar results. 

The results of our photometric comparison are shown 
in the left panels of Figure [BJ Overall, field brown dwarfs 
produce poor fits to the HR 8799 b photometry. This 
may be a result of comparing objects with an order of 
magnitude difference in surface gravity (log g~4.0-4.5 
[cgs] for HR 8799 b vs. log g~5.0-5.5 for field brown 
dwarfs). The differences can also arise from comparing 
objects with different cloud properties and/or metallici- 
ties. T dwarfs produce the worst fits as a result of their 
blue J-K colors compared to HR 8799 b. The best- 
fitting spectral type (as measured by the average G'-value 
for a given spectral type) is L5, although similarly good 
G values are obtained from L5 to L8. However, typical 
L5 objects with G values near the mean for that spectral 
type do not produce a strong match because of the red 
colors of HR 87 99 b compared to other field objects (see 
also Figure 6 of I Alters et all 120 101) . 

The best fitting individual objects produce much bet- 
ter matches to the HR 8799 b NIR photometry. The 
left panel of Figure [7] shows the five best-fitting ul- 
tracool dwarfs, which have optical spectral types be- 
tween L6 and L8 (for those with optical classifica- 
tions). These objects all have very red NIR colors 
(J-Ks~ 2. 0-2. 5 mag) and generally show evidence of 
low gravity and/or thick clouds. For example, atmo- 
spheric model fitting to the 0.8-14.5 fxm spectrum of 
SDSS J085758.45+570851.4 indicates that it has a low 
surface gravity (log q=4.5) and a very cloudy atmo- 
sphere ((Stephens et al.l 12009( 1. HR 8799 b is well fit 
by 2MASS J22443167+2043433 (2MASS J2244+2043), 
which is the reddest L dwarf known w ith a J-K$ color 
of 2.48±0.15 mag l|Dahn et all 120021) . A comparison 
of atmospheric models to mid-i nfrared photometry o f 
2MASS J2244+2043 presented bv lLeggett eFaXl (|2007bD 
suggests that it has extreme cloud properties with strong 
vertical mixing (AT 2Z >10 4 cm 2 s _1 ). There is also evi- 
dence that 2MASS J2244+2043 has a low surface gravity 
based on weak K I lines and FeH bands (|McLean et all 
120031: Hooper et all 120081) . The overall best- fitting ob- 
ject to HR 8799 b is 2MASS J21481628+4003593, which 
has a J~K$ co lor of 2.38±0 . 06 m a g and a triangula r 
iJ-band s hape (lLooper est all [200l 1 Alters et all 120071 ). 
However, lLooper et all (|2008l ) interpret this object as 
old and metal-rich. Note that although the NIR SED of 
HR 8799 b is well-matched by these red field L dwarfs, 
the luminosity of HR 8799 b is roughly 3-10 times lower 
than these best-fitting objects. 

4.3. 1.1-4-1 fJ-m Photometry 

We now extend the empirical comparison to 4 /jm by 
incorporating L'-band photometry. Wc use the compila- 
tion of u ltracool dwarf near- and mid-infrared photome- 
try from lLeggett et al.l (|2010l ) to fit the J, H, K s , and V 
photometry of HR 8799 b. A'mko magnitudes are con- 
verted to Ks magnitudes using a polynomial fit to dif- 



ferenced synthetic photometry of L and T dwarfs from 
the SpeX Spectral Prism Library (see Appendix Un- 
certainties in the polynomial fit are accounted for in the 
As-band photometric error. 

To incorporate the photometric measurement errors of 
the comparison objects into our fitting analysis we define 
a goodness-of-fit statistic G" as follows: 



(6) 



where aj^ k . is the photometric measurement uncertainty 
of the ultracool dwarf k for filter t[|3 The scaling factor 
C k is calculated iteratively as described in Appendix [Bj 
The weights w\ are defined as in Equation [5] and corre- 
spond to values of 0.11, 0.19, 0.22, and 0.48 for the J, 

H, Ks, and V filters. 

Mid- to late-L dwarfs provide the best fits to the 1- 
4 /xm photometry of HR 8799 b (Figure [H right pan- 
els). These results hold when the weights are set to 
unity for all bands. Similar to our empirical comparison 
of the NIR SED (1.1-2.4 pan), the best-fitting field ob- 
jects to the near- and mid-IR photometry (1.1-4.1 /zm) 
are late-L dwarfs (L7.5-L8 optical spectral types) with 
red NIR colors, most of which exhibit evidence of low 
gravities and/or abnormally dusty photospheres (Fig- 
ure H right panel). 2MASS J08251968+2115521 (here- 
after 2MASS J0825+2115) has the lowest G" value and 
is a well-studied red L dwa rf (</-A,s'=2.07±0.04 mag; 
L7.5 optical spectral type). iStephens et all (|2009f ) and 
iCushing et all (|2008f ) performed atmospheric model fits 
to its near- and mid-infrared spectrum and found best- 
fitting effective temperatures of 1200 K and 1400 K, re- 
spectively, with unusually thick clouds compared to other 
field L dwarfs. 

The location of HR 8799 b in the ( J-H, K s -L') color- 
color diagram is shown in Figure [5J HR 8799 b has a 
very red J-H color (1.43 mag) compared to field L and 
T dwarfs. Its J-H color is even more extreme than the 
reddest field L dwarfs (labelled in Figure [5J , which are 
generally thought to have low gravities and / or unusually 
dusty photospheres. 

In Figure [5] we compare the 1.1-4.1 (atxi photome- 
try of HR 8799 b to four well-studied, young, low- 
mass L dwarfs. AB Pic b is a ~13-14 Mj up com- 
panion to a young (^30 Myr) K 2-type member of th e 
Tucana-Ho rologium associ at ion (Chauvin et all 120051 ). 
Recently, iBonnefov et all (|2010 ) obtained medium- 
resolution integral-field spectroscopy of AB Pic b from 

I. 1-2.5 /im and determined a spectral type of LO- 
LL A comparison of HR 8799 b to the Bon- 
nefoy et al. NIR spectrum in Figure [9] shows a 
decent match from 1.1-2.4 /im, but HR 8799 b is 
slightly redder in J-A5. Similarly, HR 8799 b has 
redder NIR colors compared to the integrated-light 
spectrum of SDSS J224953.47+004404.6AB (hereafter 
SDSS J2249+0044AB), a pair of low-gravity L d warfs 
(L3+L5) recently discovered by I Alters et all ((20101) . The 
NIR spectrum of the young ~8 Mj up L4l2-type com- 

11 We remove angle brackets from J 7 /, i m Equation [6] to dis- 
tinguish between actual photometry and synthesized photometry 
from spectra, which we use in Equation [3] 
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Fig. 7. — Left: Best-fitting objects to J-, CH4S-, H-, and Xg-band photometry of HR 8799 b (at 10 pc; see also Figure Qj). These 
ultracool dwarfs generally have late-L optical spectral types, red colors compared to normal field dwarfs with the same spectral types, 
and indications of low surface gravitites and/or cloudy atmospheres. Although the NIR SEDs provide good matches, the luminosi- 
ties of these field objects are ~3-10 times higher than th at of HR 8799 b. T he SpeX/prism spectra of 2MASS J21481628+4003593 
and 2MASS J2244316 +204343 were originally published in [Cooper ct al. (200S|), SDSS J010752.33+004156.1 and SDSS J0857+5708 in 
IBurgasser et ail 1120101 '). and SDSS J161731. 65+401859. 7 in lChiu et"aH (120061 ). Right: Best-fitting objects to the 1-4 /im photometry of 
HR 8799 b. See §I]and Figure [6] for details about the fitting procedure. 



anion 1RX S J160929. 1-210524 b from lLafreniere et al.l 
2001 and lLaireniere et~aTl (pOlOf ) is significantly bluer 
than HR 8799 b, although the Kg-L' color is nearly 
identical. Finally, HR 8799 b appears to be slightly 
bluer compared to the NIR spectrum of the ~5-8 Mj up 
mid /late L-ty p e obj ect 2MASS J1207-3932 b from 
IPatience et al.l (|2010|) . None of these four low- gravity 
L dwarfs provide good templates to the 1.1-4.1 (xm pho- 
tometry of HR 8799 b. 

Altogether, our spectral and photometric comparisons 
to field brown dwarfs suggest a spectral type between L5 
and T2 for HR 8799 b. Although peculiar compared to 
most L and T dwarfs in the field, the planet's photome- 
try is consistent with the reddest field L dwarfs. These 
results imply that HR 8799 b is the lowest-mass L/T 
transition object currently known. 



5. COMPARISON TO ATMOSPHERIC MODELS 

To derive the physical properties of HR 8799 b we com- 
pare our 2.12-2.23 /im spectrum and 1.1-4.1 /jm pub- 
lished photometry to atmospheric models. We use three 
published grids of low-temperature atmospheric mod- 
els in our analysis. The models of iHubenv fc Burrowsl 
(|2007l hereafter HB07) were developed to assess the ef- 
fects of nonequilibrium chemistry on the emergent spec- 
tra of brown dwarfs and giant planets. In these objects, 
departure from local chemical equilibrium (LCE) can 
arise when vertical mixing dredges up molecules from 
warmer, deeper atmospheric layers on timescales shorter 
than chemical reaction timescales and primarily affects 
CO/CH4 and N2/NH3 chemical abundances. The HB07 
grid includes a baseline LCE model as well as three vari- 
ations of the mixing timescale, which is parametrized 
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Fig. 8. — Color c olor diagram for field L and T dwarfs from 
ILeggett et al. (2010|) compared to the HR 8799 planets. L dwarfs 
are plotted as red (early-L) to orange (late-L) and T dwarfs are 
plotted as light blue (early-T) to dark blue (late-T). The HR 8799 
planets are shown as yellow stars and objects with red J-H colors 
are labelled. HR 8799 b is redder in J-H than field L and T 
dwarfs (the large uncertainty in J-H for planet d is consistent 
with late-L dwarfs). HR 8799 c shares similar colors to field objects 
near t he L/T transition. In addition to objects from lLeggett et al.l 
H2010I ). we also overplot the position of the recentl y discovered 
youn g L dwarf binary SDSS J224953.47+004404.6AB filers et al.l 
I2010I : it-band magnitudes are converted to Kg using the relation 
in Appendix [A) . 

by the eddy diffusion coefficient K zz (with modei val- 
ues of 0, 10 2 , 10 4 , and 10 6 cm s^ 1 ). HB07 compute 
models for three CO/CH4 chemical reaction timescales 
("slow," "fasti" and ll fast2 v ). We consider the slow 
and fast2 chemical timescale prescriptions (fasti and 
fast2 are similar) spanning effective temperatures of 7 00- 
1800 K (AT eff =100 K; see lHubenv fc Burrows! [20071 for 
details) . The models are for solar metallicity abundances 
and include clear atmospheres and cloudy atmospheres 
containing 100 /im forsterite particles (Figure fT0|) . The 
surface gravity of HR 8799 b is expected to be ^4.0- 
4.7 de x given its luminos ity and an age range of 30-160 
Myrs ([Baraffe et al .1120031) . so we use the log g=4.5 mod- 
els here. 

Metallicity plays an important role in shaping the 
near-infrared spectra of ultracool objects by impact- 
ing the st rength of col l ision-induced absorp tion by H2 
(CIA H 2 ; ILinskvl [1969L IBorvsow et aLl[l997t) . CIA H 2 
is strongest in the 1.3-2.5 /urn spectral region and can 
affect the shape and amplitude of the isT-band emer- 
gent spectrum; metallicity variations are therefore di- 
rectly relevant to our OSIRIS spectrum. We use the 
grid of low-temperat ure LCE atmospheric models from 
lBurrows~et al. (2006, hereafter BSH06), which sample an 
order of magnitude variation in metallicity (i?=-0.5, 0.0, 
+0.5) for effective temperatures of 700-1800 K (Figure 
[TO]) . We consider log <?=4.5 models with clear and cloudy 
variants for our analysis. As emphasized in BSH06, we 
note that there remain many uncertainties in the meteo- 
rological physics behind cloudy models in general. This 
particular set is only a representative version of particle 
size and cloud prescription which works well for modeling 
L dwarfs. 

The red colors of HR 8799 b suggest that photosphcric 



Fig. 9.— 1.1-4.1 (im photometry of HR 8799 b (black circles) 
compared to four well-studied, low-mass, low-gravity L dwarfs 
(red). The spectra are normalized between 2.0-2.4 fim and are 
offset by a constant, and the HR 8799 b SED is normalized to 
the K s band. The NIR (~l.l-2.5 fim) sp ectra of AB Pic b 
HBonnefov et al.|[20TOl'l . S DSS J2249+0044AB ifX llcrs et aH 120101, 
and 1RXS J1609-2105 IILafreniere et al.l [20081 ; ILafreniere et al.l 
|20T(T) are bluer than HR 8799 b, although the NIR SED of 
2MASS J 1207-3932 b is r edder. L'-band ph otometry is included 
for SDSS J2249+0044 AB l|Allers et alJl201Cft . 1R XS J1609-2105b 
IILafreniere et al.ll20Toh , and 2MASS J1207-3932 b l IMohantv et al.l 
1200711 . None of these low-gravity L dwarfs provide good matches 
to the SED of HR 8799 b. 

dust may play an important role in shaping its emer- 
gent spectrum. The HB07 and BSH06 models we con- 
sider here include a single cloud prescription with mod- 
est dust content and large grain sizes. To compare 
with a more extreme case of dust formation we also use 
the solar me t allicit y Ames-Dusty atmospheric models of 
lAllard etTa l. (2001), which considers the limiting case of 
dust formation with no gravitational settling and is com- 
puted in LCE. Thirty types of spherical grains are in- 
cluded in the models with an interstellar size distribution 
from 0.00625-0.24 /im. The grid contains effective tem- 
peratures from 500-2000 K (AT cff =100 K) and surface 
gravities from 3.5-6.0 dex (Alog g=0.5). The effect of at- 
mospheric dust is to increase the gas temperature in the 
outer photosp heric layers compared to the dust-free case 
(|Allard et al][200l . This heating can result in weaker 
molecular features (like CH4 and H2O) and, in combina- 
tion with the depletion of metals from the gas and the 
higher grain opacities at shorter wavelengths, tends to 
smooth out the SED to more closely resemble that of a 
blackbody. This phenomenon is shown in Figure [TTJ in a 
comparison between the HB07 clear (red), HB07 cloudy 
(green), and Ames-Dusty models (blue), which progres- 
sively shows the impact of higher levels of photosphcric 
dust at different effective temperatures. (Note that the 
LCE HB07 cloudy models end at 1000 K.) 

5.1. Physical Properties from Spectroscopic Analysis 

We use the same \ 2 fitting procedure as we did for field 
brown dwarfs in <2] to compare the model atmosphere 
spectra to our OSIRIS spectrum. Our OSIRIS spectrum 
is not flux calibrated so for our analysis of the spectrum 
we cannot compute a radius as we do in i j5.2l using pho- 
tometry. For the same reason we do not include the 
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Fig. 10. — Atmospheric models from HB07 and BSH06 for changing T e g, log g, eddy diffusion coefficient (K zz ), and metallicity (Z) for 
clear and cloudy atmospheres. Clouds (thin dashed lines) produce dramatic effects on emergent spectra by suppressing the J-band flux 
and enhancing the A"-band flux. The shape and amplitude of the 2.2 fim region can be influenced by all four parameters. 



spectrum into a joint analysis with the photometry. The 
X 2 values for the HB07 and BSH06 models continue to 
decline with increasing T e g until they reach a minimum 
at -1300-1600 K (Figure[T21) . The overall behavior of the 
X 2 values is a result of the strength of the 2.2 /im CH4 
feature in the models diminishing with higher tempera- 
ture as the carbon balance shifts from CH4 to higher CO 
abundances. The result is a flattening the 2.12-2.23 fim 
spectral region in the models, producing better fits to our 
HR 8799 b spectrum. In contrast, the x 2 distributions 
for the Ames-Dusty models (Figure H"3"| are relatively flat 
for effective temperatures higher than —800 K as a result 
of inhibited methane formation. 

For effective temperatures between —800-1300 K, the 
HB07 slow models with clouds and high K zz values pro- 
duce better fits than clear LCE models. The x 2 distri- 
butions for cloudy HB07 models result minima —100 K 
lower than for the clear models. For the BSH06 mod- 
els, the cloudy versions tend to fit the data slightly bet- 
ter than clear models. Higher gravity Ames-Dusty mod- 
els tend to fit the spectrum better than lower gravities 
for T c g>900 K, although the actual goodness-of-fit for 
the Ames-Dusty models is worse than for the HB07 and 
BSH06 models as is evident from the x 2 values. 

The dominant opacity sources (CH4, CIA H2, and 
H2O) are generally well understood in the narrow spec- 
tral range of the observations, so systematic errors in the 
if-band region of the models are probably small com- 
pared t o the entire spectr al e nergy distribution. (See 
§3.2 of iBowler et all [20091 and IDupuv et all [2010l for a 
discussion of this problem.) We therefore use the Ax 2 
method to assign confidence limits as we did for field 



brown dwarfs in §3J except here we employ a x 2 value 
of 37.8 as a cutoff value between consistent and incon- 
sistent models (corresponding to a 3 it confidence level) 
for the HB07 and BSH6 grids because the minimum x 2 
value is 12.5. For the Ames-Dusty grid (which we con- 
sider distinct enough to merit its own cutoff value), the 
equivalent X 2 value is 54.2 (x 2 nin =28.9). If HR 8799 b 
has parameters that fall within these grids, and if the 
atmospheric models are correct, the HB07 and BSH06 
models constrain the effective temperature of HR 8799 b 
to >1100 K, and the Ames-Dusty models constrain the 
effective temperature to >800 K. 

In Figure HH we show examples of cloudy HB07 slow 
models with log K zz =6 (left), cloudy BSH06 models with 
Z=+0.5 (middle), and Ames-Dusty models (right) for 
log 5 =4.5 and T cff =900, 1200, and 1400 K. The 900 K 
models correspond to effective temperature s derived from 
evolutionary models dMarois et al.l l20"08bD but produce 
poor fits to the OSIRIS spectrum. 

5.2. Physical Properties from Photometric Analysis 

To further constrain the physical properties of 
HR 8799 b, we fit the same three sets of atmo- 
spheric mode l s to pu blished photometry. We use the 
iMarois et all (|2008b[) J, H, CH 4 s, K s , and L' pho- 
tometry in the fits. We exclude the published CH4I 
photometry because the 1.6 /jm CH4 band is too weak 
in atmospheri c models as a result of incomplete CH 4 
line lists fe.g.. lSaumon et al.|[2007llLeggett et al.ll2007aft . 
We choose not use other published photometry near H- 
band to avoid overweighting that portion of the spec- 
trum in the fitting (F160W: lLafreniere et al.ll2009D . As 
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Fig. 11.— HB07 clear (red), HB07 cloudy (green), and Ames-Dusty (blue) synthetic spectra for log g=4.5 and T e ff=1700, 1500, 1300, 
1100, 900, and 700 K. Increasing the photospheric dust content suppresses flux at shorter wavelengths, which then reemerges at longer 
wavelengths. At low effective temperatures the HB07 cloudy variant is similar to the clear version, but the extreme Ames-Dusty case is 
dramatically different with virtually no flux emitted in the J and H bands. Note that the HB07 LCE cloudy models end at 1000 K and so 
are not plotted in the 900 K and 700 K panels. 



de scribed below, 3 a upper limits at 3.3 (im and M-band 
by lHinz et al.1 (|2010f ) are also incorporated into our anal- 
ysis. Our approach is to fit the observed SED using two 
techniques to scale the predicted flux densities from the 
atmospheric models. 

5.2.1. Radius as a Free Parameter 

The first method uses the G k statistic (Equation [3J 
as a measure of the goodness-of-fit and scales the model 
values using C k from Equation SJ The model spectra 
represent the emergent flux at the object's surface, so 
the best-fitting C k value is also equal to (R/d) 2 , where 
R is the object's radius and d is its distance. Throughout 
this work we use HR 8799 b flux densities at 10 pc. R 
is therefore allowed to vary in a way that may not sat- 
isfy L=4:TrR 2 aT^ f{ . The model monochromatic flux den- 
sities for each filter arc computed the same way as in Sj4j 
After we compute goodness-of-fit statistics w e then con- 
sider the upper limits at 3.3 /jm and M-band (|Hinz et al.l 
12010ft to infer a range of best-fitting effective tempera- 
tures. The Hinz et al. upper limits are incorporated 
by considering two extreme cases. At one extreme, any 
models with flux above the upper limits are ruled out. In 
the weakest scheme, the upper limits are ignored and the 
best-fitting model is the one that provides the best fit to 
the detections only. As shown below, the two extremes 
produce very similar answers, so we adopt the range of 
models from both cases as the best fits. We note that 
between the 3.3 (im and A/-band upper limits, the upper 
limit at 3.3 /im provides the stronger constraint on the 



models and generally disfavors the cloudy versions with 
effective temperatures between 1400-1800 K. 

Cloudy models produce significantly better fits to the 
observed SEDs for both the HB07 and BSH06 models 
( Figure HI)]) . HB07 models with larger values of K zz yield 
better fits for a given effective temperature. The best- 
fitting models for the HB07 slow grid are the cloudy vari- 
ants with T e ff=1300-1500 K and log K zz =6, and for the 
fast2 grid are the cloudy models with T e fj=1400-1600 K 
and log K zz =6 (FigureHU). Based on the scaling factors, 
the radii inferred for the best-fitting slow models are 
0.38-0.50 i?j ur E3 and fast2 models are 0.35-0.44 i? Jup . 
For a given effective temperature, BSH06 models with 
higher mctallicities and clouds result in better fits. The 
best-fitting BSH06 model is the cloudy version with 
T e ff=1300 K and a corresponding radius of 0.48 i?j up 
(in this case the upper limits do not exclude the best- 
fitting model to the detections only). The best-fitting 
Ames-Dusty model to the detections has T c ff=1700 K, 
log .9=5.5, and i?=0.32 i? Jup (Figures [TO] and [Tp . All 
the Ames-Dusty models, however, do not agree with the 
upper limits. The M-band upper limit is inconsistent 
with T ff <900 K and the 3.3 /im upper limit is inconsis- 
tent with models >800 K. Overall, these temperatures 
and radii are hotter and smaller than those predicted by 
evolutionary models (800-900 K and 1.1-1.3 i?j up ) given 
the age and luminosity of the planet (see ^ Ej 

12 The standar d value for Jupiter's radius is 71,492 km 
IILindal c t al. 1981[), which is the equatorial radius at 1 bar. 

13 We also calculate luminosities for each model by integrating 
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Fig. 12. — Results from fitting atmospheric models to our X-band 
spectrum of HR 8799 b. The \ 2 values using the HB07 models are 
shown in the top and middle panels for the slow and fast2 chem- 
ical reaction timeseales. The bottom panel displays the results 
using the BSH06 equilibrium models for three metallicities. Mod- 
els with clear atmospheres are shown with solid lines while those 
with cloudy atmospheres have dashed lines. At low temperatures 
(<1200 K) the x 2 values increase as a result of stronger 2.2 (im 
CH4 absorption in the models. We consider models with \ 2 values 
<37.8 (99.73% c onfid ence level) as being consistent with the data 
(dotted line; see £15.11 )- 

In Figure [TS] wc show the 900 K atmospheric mod- 
els compared to the HR 8799 photometry, which corre- 
sponds to the temperature predicted by the evolution- 
ary models. The J-band fluxes are overprcdicted and 
the if-band fluxes are undcrpredicted in the HB07 and 
BSH06 cloudy models, indicating that even cloudier or 
more metal-rich models are needed to reproduce the ob- 
served flux. The 900 K Ames-Dusty model is a poor 
fit to the data, indicating that this limiting case of dust 
content is too extreme for HR 8799 b. 

There appears to be a systematic underestimate of 
radii inferred using this fitting technique compared to 
the radius derived using L=ATrR 2 aT^ s (Figure [19]) . The 
magnitude of the offset increases at lower effective tem- 
peratures and may result from inadequate atmospheric 
models. This suggests that caution must be used when 

the scaled synthetic spectra. We derive a luminosity log L\> \/Lq 
of — 5.2±0.1 dex, which is in excellent agreement with the Marois et 
al. value of-5.1±0.1 dex. We note that for lower effective temper- 
atures (larger radii), the luminosities are closer to -5.3 dex, while 
for higher effective temperatures (smaller radii), the luminosities 
are near —5.1 dex. 
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Fig. 13. — Results of fitting the Ames-Dusty models to our K- 
band spectrum. The greenhouse heating from the dust inhibits 
the formation of CH4, so the 2.2 fira feature is absent in these 
models and the \ 2 values remain similar over a wide range of tem- 
peratures ("-T000-2000 K). The inset displays the same data as a 
contour plot with contour levels of {1.1,1.2,1.3,1.4,1.5,1.8,2.} times 
the minimum \ 2 value (28.9 at T eff =2000 K, log 3=6.0). The dot- 
ted line shows the \ 2 cutoff value for being consistent with the 
data at the 3 a confidence level. The Ames-Dusty models imply 
an effective temperature >800 K for fits to our OSIRIS spectrum. 



inferring distances or radii from Cfc at very low temper- 
atures. 

To study the influence of surface gravity on our results, 
we also fit the HB07 and BSH06 models at log g—5.5 and 
5.0 (Figure [20]) . The distributions of G' values shift to 
lower effective temperatures at lower gravities. Addi- 
tionally, the minima of the distributions occur at lower 
G' values for lower gravities, indicating that low gravi- 
ties produce better fits than high gravities. We note that 
for the HB07 models the slow models result in lower G' 
values than fast2 models for all gravities, and for the 
BSH06 models the high metallicities are preferred over 
solar-metallicity and metal-poor models for all gravities. 

5.2.2. Fixed Radius 

Our second fitting method uses L= 4irR 2 <TT* s , the 
measured luminosity (log Lboi/ £©=-5.1 dex), and the 
model effective temperature to derive the scaling factor 
for each model. All models with the same effective tem- 
perature will have the same radius for a given luminos- 
ity, and therefore the same scaling factor Ck~{R/d) 2 . To 
incorporate the uncertainty in the measured luminosity 
(0.1 dex) we fit the atmospheric models in a Monte Carlo 
fashion by randomly drawing a new luminosity value, 
computing a new scaling factor using the model effective 
temperature, and calculating the G k value for each trial. 
We use the same weights Wi defined in Equation [H Wc 
repeat this for 10 3 trails, resulting in a distribution of G k 
values for each model. 

The outcome of fitting the models by fixing R at each 
T e ff are similar to those from treating R as a free param- 
eter (Figure 121]) . Cloudy models, high K zz models, and 
high mctallicity models produce the best fits to the pho- 
tometry. However, the systematic differences between 
the models and the data are similar to those noted in 
^5.2.11 with the largest discrepancy in the J-band. 
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Fig. 14. — Examples of atmospheric model fits to our iC-band spectrum of HR 8799 b (black circles with error bars) 
shows the effect of varying T e fj in the cloudy version of the slow HB07 models for log g=4.5, Z=0, and K zz =10 6 cm 2 
temperatures are shown in the middle panel for the cloudy version of the BSH06 models at high mctallicities (Z=+0.5) and log g=4.5. 
Effective temperatures below ~1200 K arc inconsistent with our OSIRIS spectrum for the physical parameters sampled in the HB07 and 
BSH06 model grids. The right panel shows the same temperatures for the Ames-Dusty models, which arc generally featureless in this 
spectral region. All models are scaled to the OSIRIS data using Equation [2] 
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Fig. 15. — Atmospheric model fitting to HR 8799 b pho tometry. 
Here the models are scaled by allowing R to vary f i|5,2.1l l. Models 
with clear atmospheres are shown with solid lines while those with 
cloudy at mos pheres have dashed lines, and colors are the same as 
in Figure H"2l Models inconsist ent with the upper limits at 3.3 (im 
and M-band ijHinz et alj20Tufl are marked with crosses. The best- 
fitting models have cloudy atmospheres with effective temperatures 
of 1300-1600 K, high-/\z Z values (top, middle), and high mctallic- 
ities (bottom). 
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Fig. 16.— Best-fitting HB07 slow (top), BSH07 (middle), 
and Ames-Dusty (bottom) atmospheric models to photometry of 
HR 8799 b (at 10 pc; black circles with error bars). The best- 
fitting models to the detections for the HB07 slow and Ames-Dusty 
models are incon sistent with the 3 a upper limit at 3.3 ^m from 
IHinz et al.l 1 120101 . horizontal lines) . The effective temperatures and 
radii inferred from these fits arc inconsistent with the evolutionary 
model-inferred values of 800-900 K and 1.1-1.3 -Rj up . Although 
cloudy HB07 and BSH06 models provide better fits than the clear 
models, the J-band flux is still overestimated in the models. This 
suggests that even higher dust opacities are required, although not 
as high as in the Ames-Dusty models, which are inconsistent with 
the upper limits. All models have been smoothed for bett er ren - 
dering and the radius is allowed to float in the scaling f i|5,2.1ll . 



6. DISCUSSION 

Atmospheric model fitting to our OSIRIS spectrum 
and published photometry suggests that HR 8799 b is 
an object with extreme physical parameters (see Table [2] 
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Fig. 17. — Ames-Dusty atmospheric model fitting to the photom- 
etry of HR 8799 b. The global minimum is located at log g=5.5 
and T e fj=1700 K, and a local minimum exists at log g=3.5 and 
T G ff=1600 K. The inset shows the same data as a contour plot for 
contour levels of {1.3,2,4,8,15,20,30} times the minimum G' value 
(0.69). Although the best-fitting models fit the detections rather 
well (indicated by the low G' values), all m odels are inconsistent 
with the upper limits from Hinz ct al. 2010. The Af-band upper 
limit is inconsistent with T c fj<900 K while the 3.3 fim upper limit 
is inconsistent with models > 800 K. The best-fitting model to the 
detections is shown in Figure H"6l 



Fig. 19. — Radii from fitting atmospheric models to photome- 
try of HR 8799 b. Values are computed from the scaling factor 
C'=(R/d) 2 for the HB07 models (red) and BSH06 models (gray). 
Cloud-free models are represented with filled symbols and cloudy 
models with open symbols. The radii are systematically offset from 
the values expected based on L=47ri? 2 crT*j using the luminosity 
of HR 8799 b (log L\, \/ Lq=— 5.1±0.1) and the effective tempera- 
ture of the atmospheric model (black curve and dashed error bar). 
The magnitude and significance of the offset for the cloudy models 
appear to increase at at lower effective temperatures. 
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Fig. 18. — Comparison of atmospheric models with effective tem- 
peratures predicted by evolutionary models (~900 K) to the pho- 
tometry of HR 8799 b. The J-band flux is overestimated and the 
A"-band flux is underestimated in the HB07 slow (top) and BSH06 
(middle) models, and in both cases strong 2.2 (im CH4 absorption 
is present, which is not seen in our A"-band spectrum. The 900 K 
Ames-Dusty model provides a very poor fit to the data. Models 
with intermediate clouds between the HB07/BSH06 cloudy ver- 
sion and the Ames-Dusty version may provide better fits to the 
photometry while still inhibiting strong methane absorption and 
remaining below the 3.3 fira and M-band upper limits. All models 
have been smoothed for bette r rendering and the radius is allowed 
to float in the scaling f i|5.2.1H . 
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Fig. 20. — Best-fitting model atmospheres as a function of the 
assumed surface gravity. Each point shows the best-fitting model 
to the J-, CH4S-, H-, K s -, and L'-band photometry of HR 8799 b. 
Lower gravities result in better-fitting models (lower G' values) 
with cooler effective temperatures. The HB07 slcrm models produce 
better fits than the fast2 models (left). The BSH06 models produce 
better fits and cooler temperatures at higher mctallicities across all 
gravities (right). 

and Figure [22] for a summary of our results) . The best- 
fitting atmospheric models are consistently the cloudy 
variants, which is in agreement with prev ious work in- 
dicatin g a high photosp heric dust content (|Marois et al.l 
I2008bt iLafreniere et al.ll2009t) . The unusually red color 
of HR 8799 b ( J-.ftfs=2-25 mag) is particularly striking 
given its intrinsic faintness (Mj=16.30 mag), which is 
comparable to mid- T dwarfs in the fie ld (which have 
J-K<0 mag; e.g. iKnapp et al.l I2004D. Its p o sition 
on the (Mj, J-K) diagram from iBurrows et al.l (|2006l 
their Figure 3) is near the "Case A" and "Case B" 
model tracks, which formulate clouds in a similarly ex- 
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Fig. 21.— Best-fitting HB07 and BSH06 atmospheric models 
when scaled so that the radii are fixed for a given effective temper- 
ature (defined using the luminosity of HR 8799 b: log L\, \/ Lq = 
— 5.1±0.1). Models with the same effective temperatures have the 
same radii in this fitting scheme. The results are similar to those 
when R is allowed to vary freely. G values are computed in a 
Monte Carlo fashion by varying the scaling factor G based on the 
luminosity uncertainty f EI5.2.2F . The mean G' values from each 
distribution are plotted, and for clarity we overplot only the 1 a 
error bars for the clear and cloudy LCE models. The color coding 
is the same as in Figure Il5l 

treme fashion to the Ames-Dusty models with no grav- 
itational settling of dust. However, the Ames-Dusty at- 
mospheric models are inconsistent with the upper lim- 
its at thermal wavelengths. The more modest dust 
model used in HB07 and BSH06 models overestimate 
the J and H band fluxes and underestimate the K 
band flux (which are suppressed and enhanced, respec- 
tively, in cloudy models), implying that an interme- 
diate case between the cloudy HB07/BSH06 prescrip- 
tion and the Ames-Dusty case will provide better fits 
to the HR 8799 b photometry. Other low-mass, low- 
gravity objects such as HR 8799 d (M./=15.26 mag, J- 
A" s =2.15 mag; IMarois et al.ll2008bl) and 2MASS 1207- 
3932 b (Mj=16.38 mag, J-Z vs=3.07 mag;[ Chauvin et all 
12004 iMohantv etaD 120071 ) sit in the same region of 
the color-magnitude diagram, suggesting that young 
planetary-mass objects possess radically enhanced cloud 
properties compared to every class of brown dwarf cur- 
rently known. 
Among the cloudy models, we find that those with 
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Fig. 22. — Graphical representation of our atmospheric model 
fitting results. The HB07/BSH06 model fits (green) to our OSIRIS 
spectrum imply effective temperatures >1200 K, while fits to the 
1-4 fim photometry suggest effective temperatures between 1300 K 
and 1600 K. Ames-Dusty model fits (blue) to the spectrum provide 
poor constraints on the effective temperature (>800 K). The best- 
fitting Ames-Dusty model to the 1-4 fim photometric detections 
of HR 8799 b is T oft =1700, log g=5.5, although all models in this 
grid are inconsis t ent w ith the upper limits at 3.3 fim and M band 
from l Hinz et all ((2010J). The evolutionary models of lBaraffe et al.l 
[120031 , red) predict temperatures of 800-900 K. 

high K zz values (10 6 cm 2 s~ x ) and metal-rich compo- 
sitions (Z—+0.5) are generally preferred. These results 
are particularly intriguing because field L and T dwarfs 
do not appear to have such high levels of disequilibrium 
chemistry; values of ~10 2 -10 6 cm 2 s _1 typically produce 
the best fits to near- and mid-infr ared photometry and 
spectroscopy of ult r acool dwarfs (iLeggett et al.l l2007bl: 
Saumon et all 120071 : iStephens et al.l 120091 : iGeballe et al.l 
20091) . The gas giants in our solar system have moderate 
values of K zz in their stratospheres (~10 2 -10 5 cm 2 s~ x ; 
ISaumon et al.l [2006 j3 with substantially metal-rich at- 
mospheric compositions (metal abundances range from 
~2-4 times solar values for Jupiter, up to ~10 times for 
Saturn, a nd even higher for some species in Uranus and 
Nept une; lAtreva et alj|2003t lLoddersl l2004l: iFlasar et all 
12001 iFortnev fc Nettelmannll2009D . 

The metallicity of HR 8799 b has implications for sce- 
narios of its fo rmation. Giant planets formed through 
core accretion (jPollack et al.l 119961 : lAlibert et aTll2005af ) 
are generally predicted to have metal-rich atmospheres 
as a result of planetesimal accretion, the accretion 
of me t al-rich gas, and/ or co r e erosion (lAlibert et al.l 
2005bt IGuillot fe Huesol l2006t IMousis et al.l 120091: see 
Fortnev et al.l 120081 for a thorough summary) . There is 
currently no theoretical consensus regarding the atmo- 
spheric metallicities of gas giants formed through disk 
instability, the pro posed alternative mechanism of giant 
planet formation (ICameronl 119781: IBqssI 119971), bu t re- 
cent simulations by iHelled fc Bodenheimerl (|2010D sug- 

14 Note that K zz is much higher in the convective regions of so- 
lar system giant planetary atmospheres than the radiative zone , 
reaching values of ~10 6 -10 9 cm s~ 1 HFeglev h Loddersl 119941 : 
lAtreva et aT]|T999l ; IMoses et all [2005; Sau mon et al.l 120071) . K zz 
also varies a s a fun ction of latitude and atmospheric pressure; see 
IMoses et al. (2005) for a detailed discussion applied to Jupiter's 
atmosphere. 
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TABLE 2 

Summary of Results 



Type of Fit 



Best Fits 



i<"-band Spectroscopic Analysis 



Field Brown Dwarfs 
HB07 Models 
BSH06 Models 
Ames-Dusty Models 



SpT=T2 (earlier than T4) 
T cff >1100 K 
T eff >1100 K 
T Eff >800 K 



1-4 /jm SED Analysis 



Field Brown Dwarfs 
HB07 Models 
BSH06 Models 
Ames-Dusty Models 



SpT between ~L5-L8 
T eff ~ 1300- 1600 K, cloudy, high K z 
T eff ~1300 K, cloudy, Z=+0.5 
T eff =1700 K, log g=5. 5 a 



a This is the best-fitting model to the detections, but all 
Ames-Dusty models arc inconsis tent with th e upper limits at 
6.6 (J,m and M band from lffiiiz et al.l 120101 '). 

gest that the mctallicity of the HR 8799 planets should 
be similar to HR 8799 itself if they formed in this fashion. 
Unfortunately, the mctallicity of HR 8799 is not known. 
HR 8799 is a A Bootis st ar, so while its pho t osphere is 
metal-poor ([M/H] — 0.5; iGrav fc Ka^[l99l ISadakanel 
|200H . this may be a result of recent accretion of metal- 
poor gas and its internal me tallicity could be very dif- 
ferent (see iMova et al.ll2010bl for a detailed discussion of 
this phenomenon). Ultimately, more refined predictions 
from planet formation theory combined with a better 
determination of the internal metallicity of HR 8799 and 
the metallicities of its planets will provide a key test to 
distinguish between core accretion and disk instability as 
the formation mechanising 

While our model fitting of the SED demands a very 
cloudy atmosphere, the other physical parameters are 
less well-constrained. We caution that it is not clear 
whether our results represent actual evidence for strong 
vertical mixing and high metallicity, or whether they are 
simply an outcome of using grids of models with only 
a limited variety of cloud prescriptions (cloud-free, 100 
/xm forsterite clouds, and completely cloudy). No atmo- 
spheric model that we analyzed agrees well with the SED 
of HR 8799 b, indicating that the atmospheric models are 
incorrect or that the physical parameters of HR 8799 b 
fall outside of the available grids. The relative success 
of these atmospheric models in reproducing the spec- 
tra of field brown dwarfs (e.g.. iCushing et alTl2008f ) sug- 
gests that the reason for the disagreement between atmo- 
spheric models and the photometry is most likely caused 
by a higher dust opacity in HR 8799 b than is included 
in the models. Indeed, there are many parameters in the 
cloud models that can influence the emergent spectra, 
such as particle size and clou d model prescriptio n (see, 
e.g., iBurrows et all 120061 and iHelling et afl 120081 for de- 
tailed discussions), and thus the best- fitting vertical mix- 
ing and/or metallicity value may change with a broader 

15 Note that the luminosities of HR 8799 b, c. and d do not agree 
with the l uminosities of giant plane ts formed in the "co ld start" 
models of IMarlev etaLI 120071) and IFortnev efaLI ( 120081 ') for any 
ages or planet masses. The planets are far too luminous compared 
to the models. The agreement with the "hot-start" models may 
suggest that the HR 8799 planets were formed through d i sk in - 
stability, as has been suggested by Dodson- Robinson et al. (2009) 
based on an analysis of plausible theoretical formation mechanisms. 



set of models. 

Our results highlight a significant discrepancy between 
atmospheric and evolutionary model predictions. The 
measured luminosity of planet b (log Lboi/£©=~5.1±0.1 
dex) combined with the age of HR 8799 (-30-160 Myr) 
implies an effective temperature of —800-900 a ra- 
dius of ~1.1 6-1.24 -R.Tnn, a nd a mass of —7-10 Mj up 
based on the iBaraffe et al.l ([20031 ) evolutionary models. 
However, we find the best-fitting atmospheric models 
have effective temperatures between 1300-1700 K and 
radii between —0.3-0.5 i?j U p, significantly hotter and 
smaller than the evolutionary model predictions. This 
dis crepancy for HR 8799 b has previously been noted 
bv IMarois et al.l (l2008bl). iLafreniere et alJ (j2009[ ). and 
iHinz et al.l (|2010T ). In contrast, the radii inferred from 
fitting the spectra of late-M dwarfs (which have similar 
radii of —1 i?j up ) have been in better agree ment with the 
value s predicted by evolu tionary models (iLeggett et al.l 
[20001 iBowler et al1l200l IDupuv et al.ll2009alj and like- 
W1SC for field L and T dwarfs ([Stephens et al.ll2009f ). 



The radii inferred from atmospheric model fitting of 
HR 8799 b are likely to be incorrect because giant plan- 
ets with degenerate cores cannot reach s uch small vol- 
umes at any point i n their lifetimes (e.g.. iGuillotl 120051 : 
IFortnev et al.l r2~007| ). Moreover, evolutionary models 
are less sensitive to boundary condition s (atmospheric 
parameters) than synthetic sp ectra are (jChabrier et al.1 
2000: iSaumon fc Marlevll2008D . which implicates the at- 
mospheric models as the culprit. One probable explana- 
tion is that HR 8799 b has atmospheric properties that 
are outside of the parameter space explored by current 
models. At cool temperatures <1000 K (large radii), 
the 100 /Ltm cloud prescription from BSH06 does not sig- 
nificantly improve the fit to the photometry compared 
to cloud- free models (Figures [15] and |2"TT) . However, in- 
creasing the metallicity not only improves the fit at those 
temperatures, but it also shifts the best-fitting model 
from higher temperatures to lower temperatures. Even 
higher metallicities and dustier atmospheres (but less ex- 
treme than the Ames-Dusty case) may yield a best-fitting 
model with a temperature and radius closer to the values 
predicted by evolutionary models and with methane still 
inhibited in the -fT-band spectrum. 

An alternative solution to the discrepancy between the 
atmospheric model fitting results and predictions from 
evolutionary models is that HR 8799 b is underlumi- 
nous, perhaps as a result of obstruction by an edge on- 
disk as suggested bv IMarois et al.l ([2008bD . This scenario 
is unlikely, however, because a similar underluminosity 
is seen with HR 8799 c and d. A radius discrepancy 
has also been f ound with 2MASS 1207-3 9 32 b which 
iMohantv et al.l (|2007l ) and iPatience et~all ([20101 ) inter- 
pret as an underluminosity caused by an edge-on disk. 
However, in light of the analogous problems with the 
HR 8799 planets, we favor a common origin cither from 
fitting imperfect atmospheric models, or fitting models 
that do not cover sufficient parameter space. No un- 
derluminosity is observed for the brown dwarf primary 

16 Note that the effective temperatures predicted by evolution- 
ary models based on the luminosity of HR 8799 b are relatively 
insensitive to age since the radii of brown dwarfs and giant planets 
are nearly constant with mass and age after ~100 Myr. At 1 Gyr 
the predicted tempe rature is ~1000 K and at 10 Gyr it is ~1100 K 
l|Baraffe et al.ll2003T ) 
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2MASS 1207-3932 A, which also makes this edge-on disk 
scenario unfavorable. More extensive atmospheric model 
development and fitting to the HR 8799 planets and other 
low-mass, low-gravity objects using different cloud pre- 
scriptions, metallicitics, and eddy diffusion coefficients 
may help resolve this issue. 

Our spectroscopic and photometric analysis indicates 
that the spectral type for HR 8799 b falls in the late- 
L/early-T range, corresponding to a physical regime 
where the atmospheres of field dwarfs are progressin g 
from cloud-filled to cloudless (Burgass er et al.l l2002f ). 
The fainter absolute magnitude of HR 8799 b compared 
to field ob jects appears to be an e xtensi on of the sug- 
gestion bv 'lMetchev fc Hillenbra"ndl (|2006l ) that the L/T 
transition depends on surface gravity (or equivalcntly 
age). In their study of HD 203030B, a companion to 
a young (^130-400 Myr) field star, they found a temper- 
ature of ^100-200 K cooler for this L/T transition object 
compared to field objects of comparable spectral type us- 
ing temperatures derived from evolutionary models com- 
bined with the age estimate of HD 203030. The temper- 
atures inferred for the HR 8799 planets from evolution- 
ary models arc even cooler compared to field objects (by 
300-400 K). Precise temperature derivations for brown 
dwarf binaries with dynamical mass measurements sup- 
port this idea of a gravity-dependen t L/T transition 
( Liu et al.|[200l IDupuv et all l2009bl) . iStephens et~aT1 
(|2009l ) find additional evidence bas ed on atmospheric 
model fits to field L and T dwarfs, and lSaumon fc Marlevl 
(|2008| ) make a similar suggestion based on the location 
of Pleiades L/T transition objects in color- magnitude 
diagrams. We note that the same result is obtained 
with 2MAS S 1207-3932 b; its lumi nosity (log L hol /L Q =- 
4.72±0.14; iMohantv etaD 120071 ) combined with the 
likely age of 5-10 Myrs (constrained through member- 
ship with the TW Hydrae association) implies an effec- 
tive t«n£CTatureofj^l 000 K based on evolutionary mod- 
els (jBaraffe et a l. 2003). This temperature is well below 
the L /T transition temperatu r e for field objects (~ 1200- 
1400, iGolimowski "et~atl [200l iCushing et aT1l2008l ). yet 
th e spectral type of 2MASS 1207-3932 b is m id- to late- 
L (jChauvin et al.ll2004l IMohantv et al.ll2007| ). 

7. SUMMARY 

We obtained a 2.12-2.23 /im spectrum of HR 8799 b 
with the OSIRIS adaptive optics integral field spectro- 
graph at Keck II. We performed an empirical comparison 
of field L and T dwarfs with our OSIRIS spectrum and 
previously published 1-4 /im photometry using the SpcX 
Prism Spectral Library and the compi l ation of ultracool 
dwarf photometry from lLeggett et all (|2010ft . Addition- 
ally, we fit these data using the HB07 model atmospheres, 
which include the effects of non-equilibrium chemistry 
caused by vertical mixing; the BSH06 models, which in- 
clude non-solar metallicity chemical abundances; and the 
Ames-Dusty models, which explore the limiting case of 
extreme photospheric dust content. Our results are sum- 
marized below, in Table [5J and in Figure |2"21 

• We rule out the presence of strong methane absorp- 
tion at 2.2 /im in HR 8799 b. The best-fitting spectral 
type to our OSIRIS spectrum is T2, although types ear- 
lier than T4 are consistent with our data. 

• Empirical fits of the SpeX Prism Library to the J- 
, CH4S-. H-, and ivTg-band photometry of HR 8799 b 



suggest a spectral type of L5-L8. The best-fitting field 
objects are red L dwarfs, which generally show evidence 
for low gravities and/or unusually dusty atmo spheres. 

• Empirical fits of the ILeggett et al.l (|2010D compila- 
tion of ultracool dwarf photometry to HR 8799 b J-, H-, 
K$-, and L'-band photometry suggest a spectral type 
of mid-L to late-L. Late-type, red L dwarfs provide the 
best empirical fits to the near- to mid-infrared SED of 
HR 8799 b. Similarly, HR 8799 b lies near red L dwarfs 
in the J-H / Ks~L' diagram, although it has even redder 
J-H colors than typical red L dwarfs. 

• Atmospheric model fits to our if-band spectrum im- 
ply effective temperatures >1100 K. 

• Atmospheric model fits to the J, CH4S, H, Ks, 
and L' bands favor cloudy HB07/BSH06 models with 
effective temperatures of 1300-1600 K and radii of 0.3- 
0.5 -Rjup, although no models agree well with all the data. 
The best-fitting Ames-Dusty model has T c ff=1700 K, 
log 5=5.5, and i?=0.32 i?j up , but all models in this grid 
are inconsistent with the published upper limits at ei- 
ther 3.3 /jm or M band. These effective temperatures 
and radii are inconsistent with evolutionary model pre- 
dictions of 800-900 K and 1.1-1.3 i?j up based on the lu- 
minosity of HR 8799 b and the likely age of HR 8799. 
The origin of the discrepancy is unclear, but we suggest 
it is likely a result of imperfect atmospheric models or 
inadequate range of physical parameters in the models. 

• HB07/BSH06 atmospheric models with lower sur- 
face gravities, higher K zz values, and higher metallici- 
ties generally provide better fits to the photometry of 
HR 8799 b. The effective temperature of the best-fitting 
model decreases as the surface gravity is lowered and as 
the metallicity is increased. 

• Altogether, we find that HR 8799 b has a spec- 
tral type consistent with L5-T2. With an evolution- 
ary model-derived temperature of 800-900 K, HR 8799 b 
provides further evidence that the L/T transition is a 
gravity-dependent (or, equivalcntly, age-dependent) phe- 
nomenon. 

Future photometric and spectroscopic observations can 
help to further constrain the physical properties of 
HR 8799 b. Cloudy atmospheric models indicate that 
the y-band flux should be suppressed in a similar fash- 
ion as the observed J-band flux. The M-band region is 
sensitive to non-LCE chemistry; unfortunately, none of 
the HR 8799 plan ets were detected in a 2.7 hr observa- 
tion in M-band bv lHinz "eTall (pOlOt) . The lack of strong 
CH4 absorption at 2.2 /im in our observations suggests 
that the CO feature of HR 8799 should be strong. Spec- 
troscopy targeting the _ff-band methane and the 2.3 /im 
CO feature can provide even further constraints on the 
effective temperature and nonequilibrium chemistry of 
HR 8799 b. These observations, combined with spec- 
troscopy of HR 8799 c and d, will elucidate the physi- 
cal properties of this emerging class of low-mass objects 
which are characterized by low surface gravities, low lu- 
minosities, and exceptionally cloudy atmospheres. 
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Fig. 23. — Kmxq-K $ as a function of spectral type. Differential magnitudes are computed from synthetic photometry of objects in the 
SpeX Prism Spectral Library. The black curve shows the best-fitting fifth order polynomial to the data, and the rms scatter (0.0144 mag) 
is shown as dotted curves. 
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APPENDIX 
A'mko TO K s CONVERSION 

We converted .KMKO-band to i^s-band magnitudes for our sample of L and T dwarfs from iLeggett et al.l (|2010() by 
performing synthetic photometry to objects in the SpeX Prism Spectral Library and fitting a fifth order polynomial 
to the differential magnitudes (Figure [23]) . Optical types are used for L dwarfs when available, otherwise near- 
infrared spectral types are used. The best-fitting fifth-order polynomial is given by K s~KyiKO = ^2i = i CLiSpT 1 ^ 1 , where 
a={0.05538, 0.0055409, -0.0038657, 0.00083746, -7.5248x 10" 5 , 2.0778x 10~ 6 } and SpT is the numerical spectral type 
defined such that SpT=0 for L0 (e.g., L7=7, T0=10). This relation is valid from L0 to T7, and the rms of the fit is 
0.0144 mag. 

FITTING MODELS TO DATA WITH ERRORS IN BOTH QUANTITIES 

When there arc uncertainties in both the model and the data we use the goodness-of-fit statistic G" (Equation [6]) , 

which scales both the model and the model errors by a scaling factor C k . Similar to Ck (Equation [2]) and C k 

(Equation |4{ , C k is computed by minimizing the G" merit function. Taking the derivative of G" with respect to C k , 
equating the result to zero, and simplifying, we get 



i=l 



= o, 



(Bl) 



where fi is the model flux density for filter i (of n filters), Tk,i is the flux density for object k in filter i, and c/,, and 
<jp k . are the uncertainties in both quantities. 

We solve for C k in an iterative fashion for each object with weights defined by Equation [5j C k was usually within 

30% of C k when C k was computed with the "model" photometric errors ignored. Note that when <jjr k .=0, C k becomes 

C k , and additionally when iUj=l, C k becomes Ck. 
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